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Comment on "Monte Carlo Algorithms 
with Absorbing Markov Chains: Fast 
Local Algorithms for Slow Dynamics" 

In (!]] the problem of non-critical slowing of Monte-Carlo 
simulations is addressed. This phenomenon may occur for 
instance in low-temperature simulations of the Metropo- 
lis or Glauber dynamics if the system gets stuck in a local 
minimum of the energy. To get around this problem, a dif- 
ferent algorithm for simulation of the same dynamics has 
been presented in (2J. This algorithm, commonly called 
•the n-fold way, directly simulates the times the system is 
,at rest. It is most efficient if the attraction basin of the 
energy's local minima are small, that is, if the probabil- 
ity to leave a basin after one step and afterwards to stay 
outside for a long period of time is high. 

It may however happen that the attraction basins are 
of moderate size, for example if some energetically de- 
generate configurations are in the same basin. Thus it is 
reasonable to construct a generalisation of the n-fold-way 
algorithm, which simulates the sojourn times in the basins 
and the exiting configurations without altering the given 
dynamics. 

One such algorithm is proposed in jl). This algorithm, 
however, contains one element which changes the dynam- 
ics drastically and - even worse - it does not produce the 
correct stationary distribution. Let us illustrate this by 
the following simple example. Consider a random walk £ 
in {0, 1,2,..., 100} with reflecting boundaries at and 100 
and probability | to move left or right at each time step 
inside. The stationary distribution of £ is P(£ = n) = 
for 1 < n < 99 and P(£ = n) = ^ for n = 0, 100. To 
every state n we associate a basin <B„ := {n} if n < 40 
or n > 60 and 25„ := {n,n + 1} if 40 < n < 60. For 
this setting we have performed a simulation according to 
■the recipe given in [Q. The algorithm produces a single 
limiting distribution depicted in Fig. 1, which is clearly 
distinct from the correct one. It effectively introduces a 
spurious drift depending on the basins chosen. This will 
also change dynamical properties. 

We finally give the relevant probabilities from which 
a correct algorithm may readily be constructed in the 
spirit of Consider a Markov chain £ in the config- 
uration space with given one time step transition prob- 
ability P(E(i + 1) = a' | £(*) = a) = pi(a' | a). To ev- 
ery configuration a associate a set 23 a = : {ax, ■ ■ ■ , &b„} 
of b a > 1 basin configurations containing at least a =: 
<j\ itself. With the help of the b a x ^-matrix B a := 
(jPlip'k I o~i))k,l=i,...,b a the conditional distribution of the 
first exit time 

T t := min{t + r : E(i + r) ^Q5 s(t) ;r > 1} (1) 

from the basin can be written as 

F(T t = t + r | E(t) BJBT-H, (2) 
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FIG. 1: Distribution of the random walk described in the text 
using Eq. (3) (o) or Eq. (^) (O) averaged over 10 7 time steps 
starting from a uniform initial distribution. 

where we have used the 6^-vectors ft := (1,1,..., 1), 
$ := (1, 0, 0, . . . , 0) and the dagger denotes transposition. 
The distribution of the exit times is in accordance with 
that stated in Q. Now form the set <£ a := {ai, . . . , a Pa } of 
all e a exit configurations outside Q5 CT that can be reached 
from in one time step. The conditional probability to 
be in one of these configurations at given first exit time is 


p(E(T f ; 


\T t =t + r; £(t) = a) 


/3t(l-B CT )BJ-V 


(3) 


where we have introduced the e„ x & CT -matrix E a := 

(pi(a n \5-k))n=i,...,e a ;k=i,...,K and tne e^-vector {e n ) m := 
for m 7^ n, (e„)„ := 1. The algorithm proposed in 
would correspond to substitute the right-hand side of (3) 

by 


Er = >(i 


B„)Bl'8 


(4) 
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